ADP-09-001/T679 



Stout-link smearing in lattice fermion actions 



o 
o 

< 
(N 



J. B. Zhang, ^'^ Peter J. Moran,^ Patrick O. Bowman,'^ Derek B. Leinweber,^ and Anthony G. WiUiams^ 

' ZIMP and Department of Physics, Zhejiang University, Hangzhou, 310027, P. R. China 
^Special Research Center for the Subatomic Structure of Matter ( CSSM) and Department of Physics, 

University of Adelaide 5005, Australia 
"^Centre for Theoretical Chemistry and Physics and Institute of Natural Sciences, 
Massey University (Albany), Private Bag 102904, North Shore City 0745, New Zealand 

(Dated: August 26, 2009) 

The properties of the momentum space quark propagator in Landau gauge are studied for the 
overlap quark action in quenched lattice QCD. Numerical calculations are performed over four 
ensembles of gauge configurations, where three are smeared using either 1, 3, or 6 sweeps of stout- 
link smearing. We calculate the non-perturbative wave function renormalization function Z{p) and 
the non-perturbative mass function M{p) for a variety of bare quark masses. We find that the 
wave-function renormalization function is slightly sensitive to the number of stout-link smearing 
sweeps. For the mass function we find the efi'ect of the stout-link smearing algorithm to be small 
for moderate to light bare quark masses. For a heavy bare quark mass we find a strong dependence 
on the number of smearing sweeps. 
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The quark propagator is one of the fundamental com- 
ponents of Quantum Chromodynamics (QCD). Although 
it is not a physical observable, many physical quantities 
are related to it. By studying the momentum-dependent 
quark mass function in the infrared region we can gain 
valuable insights into the mechanisms of dynamical chiral 
symmetry breaking and the associated dynamical gener- 
ation of mass. At high momenta, one can also use the 
quark propagator to extract the running quark mass [H . 

Lattice QCD techniques provide an avenue for the non- 
pertubative study of the quark propagator. There have 
been several lattice studies of the momentum space quark 
propagator I, i, i, i, 1 0, i, 1 IH EH using different 
fermion actions. Finite volume effects and discretization 
errors have also been extensively explored in lattice Lan- 
dau gauge • 

The overlap fermion formalism (l^ . [Tst realizes an ex- 
act chiral symmetry on the lattice and is automatically 
0{a) improved. There are many salient features of over- 
lap fermions, which include no additive renormalizations 
to the quark masses, an index theorem linking the num- 
ber of zero-modes of the Dirac operator to the topologi- 
cal charge Q, and evading the so called "no-go theorem" 
etc., however they are rather computationally demand- 
ing. There are many suggestions on how to reduce the 
computational cost. One such proposal is the use of a 
more elaborate kernel, together with a fattening of the 
gauge links 

The idea of any UV- filtered fermion action [l^, [2l|, 
Is^ . [35| is that one will carry out the calculation on a 
smoothed copy of the actual gauge field and evaluate the 
Dirac operator on that background. This yields a new 
fermion action which differs from the old one by terms 
which are both simultaneously ultralocal and irrelevant. 
The term "UV-filtered" indicates that such an action is 
less sensitive to the UV fluctuations of the gauge back- 
ground. Sometimes, one also speaks of "fat-link" actions. 



There is a great amount of freedom available when 
generating a smoothed copy of some gauge field. One 
needs to decide on the smoothing recipe (APE [l7[, 
HYP [l^, stout-link [l^l, etc.), on the parameter 
(ck^^^, 0:^^3,0^'^ = p) and on the number of iterations, 
^iter^ In any case, with fixed {a,v}^'^'^) the filtered "fat- 
link" action is in the same universality class as the usual 
"thin- link" version . Unfortunately, if any smoothing 
process is over-applied, some important properties of the 
theory are lost. Therefore, one needs to find a balance 
between the smoothing procedure, which will accelerate 
convergence of the quark operator inversion and improve 
the localization properties, at the danger of losing im- 
port antphysics. Recently, Stephan Durr and collabo- 
rators [l^ liO, HH applied 1-3 sweeps stout-link smear- 
ing (2^ to the lattice gauge configurations and analysed 
how this affected various physical quantities. They claim 
that it is safe to use 1-3 sweeps of standard stout-link 
smearing on the gauge configurations. More recently, 6 
sweeps of stout-link smearing was used in the Science 
article exploring the hadron mass spectrum [36j . 

In this paper, we investigate the momentum space 
quark propagator on quenched gauge configurations. We 
utilise both the original lattice configurations and also 
the configurations which are produced by one, three, and 
six sweeps of standard stout-link smearing respectively. 
We compare results across all four cases, in order to ex- 
plore the effect of smearing on the quark propagator with 
different quark masses, different lattice momenta, etc. 

The massive overlap operator can be written as [23| 
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where Hj^{x,y) = j5Dj^{x,y) is the Hermitian Wilson- 
Dirac operator, e^H^) = H^/ \pH^ is the matrix sign 
function, and the dimensionless quark mass parameter /i 
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where is the bare quark mass and m^, is the Wilson 
quark mass which, in the free case, must He in the range 
< TOti, < 2. 

The bare quark propagator in coordinate space is given 

by 



where 
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When all interactions are turned off, the inverse bare 
lattice quark propagator reduces to the tree-level version, 
and in momentum space is given by 

{S^'Y\p)^^{T.C^:\pH}j +B^"\p), (5) 

where p is lattice momentum. One can calculate 3'^°'^ (p) 
directly by setting all links to unity in coordinate space, 
doing the matrix inversion, and then taking its Fourier 
transform. It is then possible to identify the appropriate 
kinematic lattice momentum q directly from the defini- 
tion 
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The form of q^(p^) is shown and its analytic form given 
in Ref. . Having identified the appropriate kinematical 
lattice momentum q, we can now define the bare lattice 
propagator as 



Z{P) 
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This ensures that the free lattice propagator is identical 
to the free continuum propagator. Due to asymptotic 
freedom the lattice propagator will also take the contin- 
uum form at large momenta. In the gauge sector, this 
type of analysis dramatically improves the gluon propa- 
gator [23, S, [2^. 

The two Lorentz invariants can then be obtained via 
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Here Z{p) is the wave- function renormalization function 
and M{p) is the mass function. The above equations 



imply that Z(j)) is directly dependent on our choice of 
momentum q, whilst M{p) is not. 

Standard stout-link smearing, using an isotropic 
smearing parameter psm, involves a simultaneous update 
of all links on the lattice. Each link is replaced by a 
smeared link U^{x) [2^ 



where 



Uf^{x) = exp(iQ^(a;)) U^{x) 



(10) 
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with 



ilf^{x) = psm '^{^ X 1 loops involving t/^(x)} . (12) 

We work on 16'^ x 32 lattices, with gauge configura- 
tions created using a tadpole improved, plaquette plus 
rectangle (Liischer-Weisz [23]) gauge action through the 
pseudo-heat-bath algorithm. The lattice spacing, a = 
0.093 fm, is determined from the static quark poten- 
tial with a string tension of ^/a = 440 MeV [2S]. The 
number of configurations to be used for each ensem- 
ble in this study is 50. The first smeared ensemble 
is created by applying one sweep of stout-link smear- 
ing to the original configurations with a smearing pa- 
rameter of p = 0.10. The second smeared ensemble is 
created using three sweeps of stout-smearing with the 
same value of a. We work in an C'(a^)-improved Lan- 
dau gauge, and fix the _gauge using a Conjugate Gradient 
Fourier Acceleration [29| algorithm with an accuracy of 
= J2 \dfiA^{x)\ < 10~^^. The improved gauge-fixing 
scheme was used to minimize gauge-fixing discretization 
errors '30|. 

Our numerical calculation begins with an evaluation of 
the inverse of -D(/i) on the unfixed gauge configurations, 
where £>(/i) is defined in Eq. ([T|). We then calculate the 
quark propagator of Eq. ([3]) for each configuration and 
rotate it to the Landau gauge by using the corresponding 
gauge transformation matrices {Gi{x)}. We then take 
the ensemble average to obtain S^^'^°{x,y). The discrete 
Fourier transformation is then applied to S^^'^^{x, y) and 
the momentum-space bare quark propagator, S^^'^''{p)^ is 
finally obtained. 

We use the mean-field improved Wilson action in the 
overlap fermion kernel. The value n — 0.19163 is used 
in the Wilson action, which provides mu,a = 1.391 for 
the Wilson regulator mass in the interacting case @ . We 
calculate the overlap quark propagator for 15 bare quark 
masses on each ensemble by using a shifted Conjugate 
Gradient solver. The bare quark mass wP is defined by 
Eq. In the calculation, we choose the mass parameter 
/i = 0.009, 0.010, 0.012, 0.014, 0.016, 0.018, 0.021, 0.024, 
0.030, 0.036, 0.045, 0.060, 0.075, 0.090, and 0.105. This 
choice of /i corresponds to bare quark masses, in physical 
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units, of mO = 53, 59, 71, 82, 94, 106, 124, 142, 177, 212, 
266, 354, 442, 531, and 620 MeV respectively. 

The partial results for the mass function M{p) and 
the wave-function renormalization function Z^^-* (p) = 
Z(C^\p) on a 16^^ x 32 lattice without any smearing in 
Landau gauge were reported in Ref . [ill ■ Here we focus 
on a comparison of the behavior of the overlap fermion 
propagator when using different numbers of stout-link 
smearing sweeps. All data is cylinder cut [11]. Statisti- 
cal uncertainties are estimated via a second-order, single- 
elimination jackknife. 

In a standard lattice simulation, one begins by tuning 
the value of the input bare quark mass m° to give the de- 
sired renormalized quark mass, which is usually realized 
through the calculation of a physical observable. How- 
ever, smearing a lattice configuration filters out the ul- 
traviolet physics and the renormalization of the mass will 
be different. To some extent, the effect is similar to that 
of an increase in the lattice spacing a. After smearing, 
the same input rr^ will therefore give a different renor- 
malized quark mass. The input bare quark mass must 
then be re-tuned in order to reproduce the same physical 
behavior as on the unsmeared configuration. 

We wish to directly study how the quark propaga- 
tor S(j)) is affected by smearing, through a calculation 
of the mass M(j)) and wave-renormalization Z(p) func- 
tions. In order to replicate the re-tuning procedure de- 
scribed above, we begin by first calculating M{p) and 
Z{p) for all values of m" listed previously, over all four 
types of configurations. We then select a value of the bare 
quark mass m° to investigate, and force the mass func- 
tions M{p) to agree at a given reference momentum, C,. 
This is achieved by interpolating M{p), for the smeared 
configurations, between neighboring values of the bare 
quark masses, in order to determine the required effec- 
tive bare quark mass. Any reasonable choice of C should 
suffice. By reasonable, we mean any point out of the 
far infrared (IR) or ultraviolet (UV) momentum regions, 
where lattice artifacts will spoil the results. 

In comparing the renormalization function, we first in- 
terpolate Z{p) to the effective bare quark mass, obtaining 
Z^^\p). We then multiplicatively renormalize Z'-^^{p) to 
= Z(C,p), subject to ZiCO = 1- 

We begin with a comparison of the functions M{p) and 
^(R) (p) for a small bare quark mass, with three choices of 
the reference momentum C = 2.0, 3.9 and 6.0 GeV. The 
interpolated mass functions for the smallest bare quark 
mass m° = 53 MeV are given in Fig. [1] We note the 
significant reduction in the statistical error, even after a 
single sweep of smearing. For all choices of C, the mass 
functions display strong agreement over all four levels of 
smearing, with the only differences occurring in the most 
infrared points. For the function Z''^^{p) the effect of 
smearing is also subtle, however the link smearing does 
introduce a minor splitting in the UV region. This split- 
ting leads to small differences in the lower momentum 
regions of Z*^^^(p) when we select ^ = 6.0 GeV. 

Next we consider a moderate bare quark mass of 



177 MeV, for which the functions M{p) and Z'^^\p) are 
shown in Fig. [2l As in the case of a small bare quark 
mass, we find that the mass function appears indepen- 
dent of the choice of reference momentum, however the 
discrepancy at the most infrared point is no longer ap- 
parent. The renormalization function displays the same 
splitting in the UV region. The effect of smearing on the 
quark propagator still appears to be relatively minor at 
this value of m° . 

Finally we consider a larger choice of the bare quark 
mass, to" = 531 MeV. A consideration of the mass func- 
tions M{p) given in Fig.[3]reveals a strong dependence on 
the choice of reference momentum We see that a choice 
of either ^ = 3.9, or 6.0 GeV leads to large discrepancies 
in both the low and moderate momentum regions. With 
a choice of C = 2.0 GeV we are able to obtain agreement 
in the low momentum region. 

The dependence of M{p) on C, indicates that the sup- 
pression of ultraviolet fluctuations by the smearing al- 
gorithm has spoiled the physics of the theory above 
~ 2 — 3 GeV, for this value of to°. These effects are 
clearly visible after just a single sweep of smearing at this 
heavy bare quark mass. We further note that in the case 
of 6 sweeps and C, = 6.0 GeV, the mass function drops 
to the bare quark mass. This is a clear indication that 
the Compton wavelength of the quark is small enough to 
reveal the void of short-distance interactions following 6 
stout-link smearing sweeps. 

The renormalization functions Z'^^ (p) for a heavy bare 
quark mass of mP = 531 MeV are also provided in 
Fig. [21 Apart from the small splitting in the UV re- 
gion, Z^^^p) still appears to be mostly unaffected by 
the smearing algorithm. In Fig. |4] we show the differ- 
ences in Z^^'^{p) between the smallest and largest bare 
quark masses, where in order to examine the UV split- 
ting we choose ( = 2.0 GeV. Figure S] shows that the 
magnitude of the splitting in Z'-^^ (p) introduced by the 
smearing algorithm is unaffected of the input bare quark 
mass. 

The stout-link smearing procedure can save a large 
amount of compute time in the calculation of hadronic 
physics. Not only is the Dirac operator easier to invert 
but statistical errors are reduced significantly. The con- 
clusion drawn from this study is that up to six sweeps of 
stout-link smearing sweeps induces rather small effects 
on the quark propagator for small and moderate bare 
quark masses, as claimed by Durr, et al. 19,^,^]. Af- 
ter an appropriate rescaling of the bare quark mass, the 
renormalized quark propagator displays the same physics 
as the untouched configuration. The only notable excep- 
tions are order 2% discrepancies in the renormalization 
function for all quark masses and the most infrared point 
of the lightest quark mass function. There an effect ap- 
proaching 2(7 is revealed. 

These subtle effects provide some evidence of a link 
between small topologically nontrivial gauge field config- 
urations linked to dynamical chiral symmetry breaking 
through their production of approximate zero-modes in 
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FIG. 1: The interpolated mass M{p) and renormalization Z^^-'(p) functions for the small bare quark mass, m° = 53 MeV, 
with three choices of C,. The effective bare quark masses are given in square brackets. There is good agreement in M(p) for all 
choices of ^ with up to six sweeps of stout-link smearing. A small splitting in the UV region of Z'^' (p) is apparent after three 
sweeps of smearing. This leads to a disagreement in Z^'^\p) for a large choice of = 6.0 GeV. 



the Dirac operator. Upon smearing this short distance 
physics is modified. 

Certainly the effects are subtle. However, they may 
require further investigation in the event that fermion ac- 
tions, in which all links of the action are smeared, become 
the action of choice for calculating the physics beyond the 
standard model. 
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FIG. 2: The interpolated mass M{p) and renormalization Z'^-'(p) functions for the moderate bare quark mass, mP = 177 MeV, 
with the three choices of (. The effective bare quark masses are given in square brackets. As with the small bare quark mass, 
the mass function displays good agreement for all choices of and there is also a small splitting apparent in the UV region of 
Z^'^^p). We note that the differences in Z^'^\p) appear to be independent of the bare quark mass. 
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FIG. 3: The interpolated mass M{p) and renormalization Z^^^ (p) functions for the heavy bare quark mass, mP = 531 MeV, 
for the three choices of ^. The effective bare quark masses are given in square brackets. We see that for this value of m", the 
choices ^ = 3.9, and 6.0 GeV lead to large differences in the moderate and infrared momentum regions of M(p). This indicates 
that the physics above approximately 3 GeV has been spoiled by the smearing algorithm. In we again find that the 

stout-hnk smearing algorithm introduces a small splitting in the infrared region. 
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